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Abstract. For any classical statistical-mechanics model with a discrete state 
space, and endowed with a dynamics satisfying detailed balance, it is possible to 
generalize the Rokhsar-Kivelson point for the quantum dimer model. That is, a 
quantum Hamiltonian can be constructed (on the same state space), such that the 
ground state wavefunction coincides with the classical equilibrium distribution. 
Furthermore the excited eigenstates correspond to classical relaxation modes, 
which (in cases with a symmetry or conserved quantity) permits extraction of 
the dispersion law of long- wavelength excitations. The mapping is natural mainly 
when the states have equal weight, as is typical of a highly frustrated model. 
Quantum and classical correlation functions are related by analytic continuation 
to the imaginary time axis. 



1. introduction 

The quantum dimer model T" of Rokhsar and Kivelson (RK) was inspired 
by the resonating- valence bond state of a quantum antiferromagnet, which had 
just been proposed as a starting point for the explanation of high-temperature 
superconductivity. That was not to be, but models of this type are of interest among 
the "highly frustrated systems" T characterized by massive degeneracy (or near- 
degeneracy). Generalizations of the RK construction were used to construct the first 
concrete lattice models in dimension d > 1 that manifestly exhibit fractionalized 
excitations HjEHlElE] 

The quantum dimer model had one nontrivial parameter V/t (see Eq. ()2.1|l . 
below); when it takes a special value ("RK point"), RK showed the exact ground 
state wavefunction is an equal-weighted superposition of all dimer coverings. This 
had the same (critical) static correlations as the classical dimer ensemble, an exactly 
solved model [3 |H| ■ Thus the RK point is a rather special kind of quantum critical 
point. 

Later the author noticed that, at the RK point, excited eigenstates correspond 
exactly to relaxation modes of the master equation for the natural Monte Carlo (MC) 
dynamics of the classical dimer ensemble |^ ll(J|. Furthermore, the dimer ensemble 
(on bipartite lattices) has a natural coarse-graining via the height representation , 
whereby it maps to a 2-1-1 dimensional interface model in its rough phase. Since the 
classical dynamics is easy to grasp, this mapping delivered the dispersion law for the 
quantum dimer model's elementary excitations at the RK point. (This dispersion was 
already understood variationally P). 

The purpose of this paper is to generalize the RK construction to any (degenerate) 
classical ensemble, showing the correspondence for several simple classical models. I 
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begin (Sec. [21) by reviewing basic notions of the quantum dimer modeim^] including 
the way that the classical and quantum models are connected, and the way we 
can comprehend the quantum dispersion if we know the classical dispersion. Any 
discrete classical model (e.g. ^21) with a dynamics satisfying detailed balance can be 
"Rokhsar-Kivelsonized" to produce a quantum model with the same mapping of the 
eigenfunctions to classical dynamics (SecOJ. 

Some examples are (i) Chakravarty's quantum 6-vertex model |13| . (ii) A classical 
Ising chain with spin-exchange (Kawasaki) dynamics, for which the RK model is the 
spin- 1/2 Heisenberg ferromagnet. (iii) A spin- 1/2 Ising model on the pyrochlore 
lattice jJB). The examples are built on a large basis set of essentially degenerate 
states. In that sense they are "highly frustrated" models ^31) whether or not the 
massive degeneracy arises from competing interactions. 

Finally (Sec.0J it is verified that the classical and quantum correlation functions 
are related exactly by a rotation of real time into imaginary time (and thus might be 
extracted more easily from simulations). 



2. Classical and quantum dimer model 

The Hilbert space of the quantum dimer model consists of all complete dimer covering 
configurations. Its Hamiltonian is customarily written 

n^-tJ2{\n) (=1 + h.c.) + y ^ (|i I) (I i| + |z) (z|) . (2.1) 

for a square lattice. In the t term, by an abuse of notation, "|l l) (z|" actually runs over 
1/3) (aj for every possible pair of configurations (/3, a), such that /S differs from a only 
by the replacement of a vertical pair by a horizontal pair of dimers on one plaquette. 
This is the elementary "flip move" for this model - the smallest possible change that 
turns one valid configuration into a different one (since the same four vertices are 
covered in either state). (For a dimer covering on a general bipartite lattice, flippable 
plaquettes are those around which every second edge has a dimer, and the flip move 
exchanges covered and uncovered edges.) The V term is diagonal in this Hilbert space, 
and can be rewritten VAffHp, where A/flip (a) (dependent on the configuration) is the 
number of flippable plaquettes in configuration a. 

RK noted that when V ^ t (the "RK point" ) , the ground state wavefunction is 

1*0) = ;^ El")' (2.2) 

where the sum is over all valid dimer configurations. Thus, the probability weight is 
Prob(a) — Afs^^, the same for each of the Afs states, just as in the classical ensemble. 
Any discrete classical model's dynamics is described by the master equation 

Mr)= E {Waf3Pf3iT)-Wpo.Pc.{T)y (2.3) 

/3[/3#a] 

Here Paij) is the instantaneous probability to be in configuration a at time r, and 
Waf3 is the transition rate to state a, given the system is in state /3. 
If we define 

I^„„ = -r„ = - J2 ^/3"' (2.4) 

f3[f3^a] 
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(total transition rate out of a), then ]2.'6\i can be rewritten as a matrix equation, 

p(t) - Wp(t). (2.5) 

Then the time evolution can obviously be decomposed into eigenmodes of the W 
matrix, labeled by eigenvalue —A where A > 0: 

Po(r)=^CAe^^^0i^) (2.6) 

A 

where (j)^^ is the (normalized) eigenvector. 

We actually are concerned only with the states that are connected by flips to 
a given initial state; within this component we can invoke a variant of the Perron- 
Frobenius theorem to assert that there is a unique steady state distribution, P!f\ 
characterized in matrix notation by WP'^*'^ = (this distribution is just t/)*-"-' with the 
sum of the components normalized to unity, not the sum of squares.) 

Furthermore, say the classical model has a Hamiltonian 'H{a) (we take a 
dimensionless TL which has already been divided by the temperature). Then the 
steady state should be the Boltzmann distribution, 

FW=e-^(")/Z, (2.7) 

where Z = e"^'^"-' is the partition function. Observe that Ti and Ti are Tioi related, 
at least not in the usual sense of taking the classical limit of a quantum dynamics. 
Of course, the dynamics should satisfy detailed balance, 

W^^P^'^ ^W^pPll'l (2.8) 

For the quantum dimer model, and also the generalized Rokhsar-Kivelsonized 
models, we will specialize to 

n^O, (2.9) 

i.e. the allowed configurations are all degenerate. In this case, P^^'>{a) is the same 
for every a and so H2.8I) reduces to saying W^a — i-e. the rate matrix is 

symmetric. Next observe that if the classical model's flip rate is always w (whenever 
a flip is possible), i.e. Wa/s — w or 0, then the matrix elements from H2.1|l are 
Tiap = —^Wai3 + ^^a^ap- Thus, at the RK point, the quantum Hamiltonian matrix 
is proportional to the classical rate matrix: 

H = --W. (2.10) 

w 

It follows, of course, that all the eigenvectors of the quantum matrix are the same as 
those of the classical matrix, and the eigenenergies are given by 

Ex = -A. (2.11) 

w 

This is the key result of the present paper. The quantum ground state eigenfunction 
H2.2|l is just the special case which has A = 0: its identity with the classical stationary 
state, noted originally by RK PP, follows since P'^*'^ is a null vector of W. 
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3. Generalizations and examples 

Now consider any classical model with discrete configurations with null Hamiltonian 
l|2.9() . Define a set of allowed "flips" (connecting two configurations) and endow the 
model with a Monte Carlo dynamics (in continuous time) such that every possible flip 
has rate w. Furthermore, on the Hilbert space {|a)}, define a quantum- mechanical 
Hamiltonian that includes those same flips (with amplitude t), as well as a term 
^■A/flip that penalizes a configuration once for every possible flip move. At the RK 
point t = V , the Hamiltonian and master-equation matrices are proportional, so as 
before the eigenvectors are the same and (|2.11|l holds. 

In many cases, one can construct a coarse-grained field from these configurations 
and infer the relaxational dynamics of the classical model at long wavelengths. Thus 
by the mapping, one also understands the low-energy excited states of the quantum 
model, a nontrivial problem if it were approached directly. 

3.1. Discrete models with "height" representations 

The quantum dimer model belongs to a class of models that are coarse-grained via a 
microscopic mapping of each microstate to an interface {/i(r)} in an abstract 2 -f 1 
dimensional space |??l llll[T7j . The interface is in its roughened phase, so the classical 
model is described by an effective free energy density oc |V/i(r)p. The normal modes 
of the classical master equation are simply capillary modes of this interface with 
eigenvalues A(q) oc |qp. Hence |H1 ^3 the quantum excitations are bosons with 
dispersion 

?i^(q) cx |q|2. (3.1) 

Height models have conserved quantities, often called "winding numbers" but 
most transparently understood as the components (V^/i, Vy/i). Local fiip moves 
cannot change the global "interface" slope V/i. Hence the configuration space is 
partitioned into subspaces, each of which has a steady state under the classical 
dynamics and a corresponding RK quantum ground state with zero energy. So, just 
as a Goldstone mode follows from a symmetry, the dispersion (|3.1|l is related to the 
degeneracy under changes of V/i. 

The earliest quantum height model {not then recognized as such) was the 
Anderson- Fazekas jlQ, approach to the s=l/2 triangular lattice antiferromagnet from 
the Ising limit. The basis is the Ising ground states, which essentially map to the 
dimer coverings of a honeycomb lattice. (The triangular Ising antiferromagnet in a 
transverse field maps directly to the honeycomb quantum dimer model A more 
interesting spin model that maps to the square-lattice quantum dimer model is the 
nearly Ising s = 1/2 antiferromagnet on the "checkerboard" lattice I18|.) 

S. Chakravarty has introduced a "quantum six-vertex" (Q6V) model ^Sli a 
generalization of a speculative state in which spontaneous orbital currents develop 
along the lattice edges, as was proposed for the pseudogap phase of cuprates. The 
currents define arrow variables, which satisfy an ice rule at vertices. The minimal flip 
move is to reverse all four arrows around a plaquette, provided they all point in the 
same clock sense. | The 6-vertex configurations have a height represensation 20 , and 
the whole T = behavior is strictly parallel to the QDM of RK 1 (assuming that 
analogous terms are included in the Hamiltonian!) In particular, when = 0, a "flat" 

X This is equivalent to the transverse field with projector in (16) and (17) of Ref. 1131 . 
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phase with gapped excitations occurs; in it, the arrows has long-range ahernating 
order of the so-called "d-density wave" type though with fluctuations. At the 
RK point V = t, however, a critical phase is predicted with dispersion Ij^.tf) . Quantum 
height models at the RK point will be further discussed in [18|. 



3.2. Spin exchange 

Say that our classical model is a set of Ising spins on a lattice of N sites in any 
dimension, with a zero Hamiltonian. Adopt the spin-conserving "Kawasaki" dynamics, 
i.e. the flip move is to exchange any nearest-neighbor pair. 

It might appear that this example is so trivial that it is not worth the observation 
that it may be considered a Rokhsar-Kivelson model. Some of the reasons why it 
may be of interest are (i) it adds weight to the conjecture that the RK point is 
typically a critical point and that the dispersion there is generically if there is a 
conserved quantity; (ii) a one-dimensional chain of this sort is the simplest example 
of the extension of the RK idea to the case that the classical states are unequally 
weighted (Sec. 13.41 below); (iii) the kagome quantum dimer model of Misguich et 
al ^1 maps to this model (see below). 

The corresponding spin- 1/2 quantum Hamiltonian is 



I ll){lU + h.c))+vJ2{\ Ti)(Ti l + UT)(iT l)-(3.2) 

(with same abuse of notation as in H2.1|l ). Converting to the notation of spin operators, 
each term becomes 

- t{StS- + 5-5+) + V[{^ + 5f )(i - 5|) + (i - 5f )(i + 5|)] , (3.3) 
where i,j are nearest neighbors; thus 

= -J^iS^S^ + SfS^) - MS^S^ - i), (3.4) 

where = 2V and J± = 2t. This is the ferromagnet with "XXZ" exchange anisotropy 
and the RK point here is the isotropic Heisenberg chain. 

Let's adapt (|2.2() to this case. The sum over states runs over every possible 
sequence of up or down spins; thus the A4 = 2^ terms can be grouped as a direct 
product of single-site terms, namely 



the ferromagnetic state with moments aligned in the +x direction. 

Since the model is isotropic, we know there are degenerate ferromagnetic states 
with moments in other directions. § This model conserves spin, so the quantum 
ground state is degenerate, as explained at the end of Sec. HJ Indeed, the sum in a 
wavefunction H2.2|l should only run over the mutually accessible conflgurations with 
(1 -|- cos6')A^/2 up-spins and (1 — cos9)N/2 down-spins. In the thermodynamic limit 
this is essentially the direct product generalizing H3.5|l . 

N 



cosi^l T)+sini0U)) , (3.6) 



§ For an alternative viewpoint on \^o)> note that 12. 8i is perfectly compatible with "H = ~h^2i 
in place of 12.91 '): then P^^\a) directly maps to 13. 6i : in general, "H could be any conserved quantity. 
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a coherent state in which aU spins are rotated from the z axis by an angle 9. 

Now, the classical relaxation dynamics is simple diffusion. The diffusion constant 
is easily obtained if we re-imagine the dynamics as exchanging all neighbor pairs 
of spins at a rate w. (If both spins point the same way, this has no effect.) For 
example, on a chain a marked spin executes a simple random walk with a total 
hopping rate of w to the right and w to the left. The long- wavelength behavior 
is diffusion, da{x,T)/dT = Dd^a{x,T)dx'^ , where (t{x,t) is the spin density, and 
D = w. The eigenvalues of this diffusion equation are X(q) = Dq^ = wq^ for small 
wavevectors q. Hence, via 12.10|l . the quantum model's excitations have dispersion 
hu>{(i) = tq^ = ^J±q^. But this is just the familiar formula for ferromagnetic spin 
waves! 

The same classical model, if endowed with a single spin flip dynamics, maps under 
Rokhsar-Kivelsonization to noninteracting spins in a transverse field - a trival, gapped 
quantum model. Now, one of the most interesting RK models is that of Ref. Its 
Hilbert space consists of all dimer coverings of the kagome lattice, which (it has been 
shown) correspond one-to-one to the possible Sz spin configurations on the triangular 
Bravais lattice (modulo a global spin reversal). Every "hop" in the Hamiltonian of 
|15| rearranges dimers around one hexagon, and this simply corresponds to flipping 
exactly one Ising spins: this model is precisely that transverse-field model. (Since 
every state is flippable, the V term is trivial in this case). 

The above insight, that (|3.2|l is the RK map of the classical ferromagnet, suggests 
a modification of the Kagome dimer model, by adopting a hop move that exchanges 
two "Ising" spins (corresponding to a dimer rearrangement around two hexagons). But 
- since the conservation of Sz has no particular meaning in the dimer geometry - we 
could make a similarity transform Hem UHfmU, where U is Sj for a particular 
y, or a product of such factors for any set of sites. Thus one has a family of 2^ 
distinct quantum dimer models, each with a hop move rearranging two hexagons, and 
each with gapless excitations labeled by a q^ dispersion; it might be interesting to 
investigate these models. 

A similar model to H3.2|l that may be Rokhsar-Kivelsonized is a classical, 
noninteracting lattice gas. In the same fashion as the above spin model, it maps 
to a quantum model of hard-core bosons with a nearest-neighbor attraction V. In fact 
it is that spin model, using a well-known correspondence in which up (resp. down) 
spins are transcribed to occupied (vacant) sites. 

3.3. Pyrochlore model 

Ref. |16| have developed a quantum spin-1/2 model on the pyrochlore lattice, which 
is the Rokhsar-Kivelsonization of the ground state ensemble of the pyrochlore Ising 
antiferromagnet. It is well known that those configurations map to those of the 
diamond-lattice ice model, with arrows along lattice edges. The appropriate order 
parameter for long wavelengths is the polarization, the coarse-graining of the ice- 
model arrow field, which is analogous to Vh in the height models and is conserved in 
the dynamics. Thus the classical dynamics is described by a diffusion equation and 
the dispersion is at the RK point. This contrasts with the oj ^ |q| dispersion 

elsewhere in the phase diagram, which Ref. have called "light." It indicates that, 
even in this model, the RK point is a quantum critical point. The q^ dispersion is 
generic at RK points, if there is a conserved quantity in the classical model, even in 
the absence of a height representation. 
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3.4- Classical dynamics with nontrivial Hamiltonian 

The RK mapping is possible even when the classical ensemble has unequal weights. 
The matrix W having elements 

W^, ^ Pi°)"'/V.,p(°)'^' ^ W'.,ei[^(")-^(«l (3.7) 

is symmetric, on account of detailed balance (|2.8|l : furthermore, since W is a similarity 
transform of W, they share the same eigenvalues. The eigenvectors are related by 

= 4>i^^ /\J Pa\ where refers to a right eigenvector of W. The corresponding 
quantum Hamiltonian matrix must be proportional to W. 

Let's work through a case of Metropohs dynamics, letting Wap = w if H(a) < 
Ti-iP), or w exp[—H{a) + H{/3)] otherwise. Hence by (I3.7|l . for /? 7^ a, 

W^„^ = u;e-ll«(")-«(«l. (3.8) 

The quantum hopping matrix element {t term) must be proportional to Wap. This 
depends only on the immediate environment of the flip location, provided that Ti. is 
a sum of local terms. However, the t term appears elaborate even in the simplest 
case, the model of Sec. 13.21 on a one-dimensional chain with H = —K'Y^SfS^j^i . 
Flips that change H are multiplied by e~'^'/^. Eq. (|3.4|l is replaced by Hym = 
-{Ja. + J±'5f_i5f+2) {SfSf+, + 5f 5f+i) -(J, + J.'^tiSf+2) (5f 5f+i - i) where 
Ji_ = <(l+e-l^l/2), J, = y(l+e-l^l/2), Ji' = t(l-e-l^l/2) and J,' = ^(l-e-l^l/^). 



4. Dynamic correlations 

An amusing (perhaps useful) corollary of Eq. I|2.1()|l is that for any generalized RK 
model, one can relate the quantum correlation function 

Cba{t) ^ (5(r)i(0)) (4.1) 

to the similar classical one C^'^'^^^ ba{t)- For the latter to make sense, the operators 
must be diagonal in Hilbert space, (a|v4|/3) — AaSafi and similarly B. 

Quite generally in a classical discrete system (starting in equilibrium) 

C^i--s^(r) = ^S^p/3(r|a)A„pW. (4.2) 

a 

where pp{T\a) is the conditional probability, given that the state at r = was a (which 
had probability Pa^). First, pp{0\a) = Sap; then (|2.6(l saysp^(T|Q;) = E CAe~'*''^(?I)^^^ 

and orthonormality of the eigenvectors (since W or W is symmetric) implies ca = ^i^'* ■ 
The latter formula (with a right eigenvector) holds even for the case of Sec. 13. 41 where 
the master-equation matrix W is nonsymmetric. 
Using Pa^ — 1/Ns, the final result is 

C^i--B^(r) =^e-^^iAPA (4.3) 

A 

where Ax = A^r'/' E„ ^a'/'a''^ (similarly Bx). 

A quantum correlation function is computed using a similarity transform by the 
time-evolution operator to convert B{t) — > e'^'^Be"''^'^, an operator acting at the 
same time as and taking the expectation in the ground state wavefunction: 

Cba{t) = (0|e'^^Pe-'^^i|0). (4.4) 
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Insert a complete set of states ^ |A)(A| on either side of the exponential between B 
and A, and note that (A|A|0) — Ax (similarly for B). This equality follows from 
|A) = (And it holds even in the case that the state weights are unequal, 

Sec. 031) We obtain 

Cba{t) ^^T. e^'^^^-^"^"iASA - C^'^^^BA (^r) (4.5) 

where I used (|2.11|l . Thus, the quantum correlations are the classical correlations in 
imaginary time. An obvious application is that any dynamic correlation, measured 
via classical Monte Carlo with sufficient precision, may be converted by analytic 
continuation to a quantum correlation function without the need to understand 
or compute the eigenfunctions. (Normally, for non-RK systems the same classical 
correlation function must be obtained from a quantum Monte Carlo simulation, based 
on a path integral and carried out in a space one dimension higher.) 

It is fairly surprising to obtain the simple correspondence (|4.5|l . It is true that the 
two evolution equations do have a corresponding time dependence, e^^'^ and e~''^^^^^ 
respectively, and the classical and quantum eigenfunctions are the same. Yet the 
quantities which evolve according to these exponentials, and which are represented 
in the vector space spanned by those eigenfunctions, are probability deviations on the 
classical side, but amplitudes (which must be squared to obtain probabilities) on the 
quantum side. 
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